Origin of Power Laws for Reactions at Metal Surfaces Mediated by Hot Electrons 
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A wide range of experiments have established that certain chemical reactions at metal surfaces 
can be driven by multiple hot electron mediated excitations of adsorbates. A high transient density 
of hot electrons is obtained by means of femtosecond laser pulses and a characteristic feature of such 
experiments is the emergence of a power law dependence of the reaction yield on the laser fluence 
Y ~ F n . We propose a model of multiple inelastic scattering by hot electrons, which reproduces 
this power law and the experimentally found exponents of several experiments. All parameters are 
calculated within Density Functional Theory and the Delta Self- Consistent Field method. With 
a simplified assumption, the power law becomes exact and we obtain a simple and very useful 
physical interpretation of the exponent n, which represents the number of adsorbate vibrational 
states participating in the reaction. 

PACS numbers: 82.53.St, 34.35.+a 



Hot-electron induced femto chemistry at surfaces 
(HEFatS) is a class of chemical reactions, where the 
energy to overcome the reaction barrier is provided by 
energetic ("hot") electrons. These electrons are often 
photoelectrons emitted from the surface when struck 
by intense laser pulses, but other possible hot-electron 
sources are Scanning Tunneling Microscope tips and 
Metal-Insulator-Metal nanodevices For these reac- 
tions, the reaction rate usually scales as a power law with 
the electron flux. For example, in a pioneering study of 
NO on Pd(lll) 0, it was shown that femtosecond laser 
pulses could induce desorption of NO, and a mechanism 
by multiple electronic excitations was identified The 
yield depended on the fluence as Y oc F n with n ~ 3.3. 
Subsequently, desorption induced by femtosecond lasers 
has been demonstrated for several other adsorbate sys- 
tems 0, 0, 0, showing non-linear yields, which can be 
fitted to power laws with exponents 3 < n < 8. It has 
also been shown that femtosecond laser pulses can induce 
surface hopping Q and oxidation reactions 0, El ] all 
of which can be characterized by power laws. 

A popular theoretical approach to the interaction of 
adsorbates with a high density of hot electrons is the 
electronic friction model [12[ where the excited electrons 
are assumed to equilibrate rapidly to a thermal distri- 
bution with electronic temperature T e . For sufficiently 
large timescales the propagation of the adsorbate is well 
approximated by semiclassical Langevin dynamics with 
an electronic friction rj e depending on T e as well as the 
adsorbate coordinates. While the friction model has cer- 
tainly been successful in reproducing various experimen- 
tal observations, it cannot account for the physical origin 
of the power law exponent n. Furthermore, measure- 
ments on the distribution of hot electrons in gold excited 
by subpicosecond laser pulses imply thermalization times 
up to 1 ps [13[, which is on the order of reaction times, 
and there are examples of laser induced surface chemistry 
where the assumption of a thermal distribution of hot 



electrons is in direct conflict with observations 0, Qj], G3 • 
In this letter we introduce a general model for first- 
principles calculations of the rates of HEFatS processes, 
regardless of the source of hot electrons. The central ob- 
ject of the model is a non-adiabatic Hamiltonian, which 
is used to calculate the hot electron induced vibrational 
transition probabilities of adsorbates. All parameters in 
the Hamiltonian are obtained from density functional 
theory and we show that the model reproduces exper- 
imentally observed power laws. Finally we will make a 
simple approximation for the transition probabilities and 
show that the power law then becomes exact in the limit 
of large electron flux and that the exponent is given by 
n = Er/Hco where Er is the energy barrier of the re- 
action and hoo is the energy quantum of the vibrational 
mode dominating the energy transfer. 

To analyze HEFatS rates involving multiple hot elec- 
trons, we consider a quadratic potential energy surface of 
an adsorbate with vibrational states \n) coupled to a lo- 
calized electronic resonance \a). The probability for a hot 
electron with energy s to inelastically scatter on the lo- 
calized state and induce a vibrational transition m — > n is 
denoted P rnn (e) and we assume a constant flux Jo of hot 
electrons incident on the adsorbate. It is further assumed 
that each vibrational quantum has a fixed lifetime T v n> 
and that there exist a maximum quantum number ur 
such that a reaction will proceed immediately if n > ur. 
The probability that one vibrational quantum survives 
the time interval At = 1/Jo between subsequent scat- 
tering events is then e ~ At / T vib t Each incoming electron 
will thus change the distribution of adsorbate vibrational 
states Q(n) until an equilibrium is reached between decay 
and reexcitation. 

If the adsorbate is initially in the ground state, the dis- 
tribution of vibrational states induced by the first elec- 
tron with energy si is 



L (n;ei) = Pon(ei). 



(1) 




The probability of the adsorbate being in the n'th vibra- 
tional state after the second electron has scattered is 

T IR — 1 

Q2(n;e 1 ,e 2 ) = ^ £i)P mn (£ 2 ), (2) 

m=0 

where pi(m) is the probability that the adsorbate is in 
the state m after the time interval At given by 



n R -l 



l—m 



Pi(m; £l ) = Y, Qi(^i)^J(e- At/T ^) r 

X (1 — e - At / T viby~'m 



(3) 



We exclude terms with m > ur since such excitations 
would have led to a reaction by assumption. Proceeding 
like this, the probability Qs(n; £i, £2, £3) of being in the 
n'th excited state after the third scattering event can be 
expressed in terms of Q2(n; £1, £2) and so forth. Since all 
vibrational states n >ur lead to a reaction, the reaction 
probability of the /c'th electron is 



tdR 

-Lb 



(4) 



For large k this will approach a limiting value, Pr. 

To calculate the vibrational transition matrix P mn (s) 
we consider a Newns- Anderson type Hamiltonian with 
substrate states \k) and a resonant state \a) linearly cou- 
pled to a number of vibrational modes with creation op- 
erators for vibrational quanta b\ [HI, [l6[ : 

H = € c{c a + e k C l C k + Yl {y^c\c k + V* k c\c^j 
k k 



FIG. 1: (Color) The difference between excited and ground 
state densities for CO adsorbed on Cu(lll). Black balls are 
Carbon atoms and red balls are Oxygen atoms. Blue contours 
is excess density in excited state and green contours is excess 
density in ground state. The excited state is constructed by 
occupying a 7r* orbital of CO. For clarity we only show the 
density difference in a single supercell. 



Self- Consistent Field designed to handle molecular or- 
bitals hybridized with metallic states and calculates the 
expectation values of excited states which are not eigen- 
states of the Hamiltonian but involve an occupied reso- 
nance. Details on the method and comparison with ex- 
periment can be found in Refs. [17| and [1 31 ] - 

As an example we calculate the fluence dependent 
transfer of energy from hot electrons to a CO molecule 
adsorbed on Cu(lll), mediated by excitation of an un- 
occupied 7r* orbital. CO adsorbs with the molecular axis 
perpendicular to the surface and the symmetry of the 
adsorbed molecule thus only allows a linear coupling to 
the center of mass (COM) and internal stretch vibra- 
tions of CO. The parameters for CO chemisorbed at a 
Cu(lll) bridge site were calculated with the code gpaw 
which is a real-space Density Functional Theory 



The model essentially describes a harmonic oscillator, 
which is displaced when the state \a) is occupied and 
the coupling V a k to the metallic states introduces a fi- 
nite lifetime of \a). As previously published [17], the 
transition probabilities P mn (£) can be calculated exactly 
in the wide band limit where the density of states pro- 
jected on the localized state \a) is a Lorentzian cen- 
tered at £0 with full width at half maximum given by 
r = 2ir^2 k \V ak \ 2 6(e () - e k ). The probabilities P mn (£) 
depend on the dimensionless parameters: huJi/Y and 
9i — (Xi/hcj) 2 and the reaction probability P R also de- 
pends on these parameters in addition to the reaction 
quantum number ur ~ Er/TujO. The quantities Er and 
uj can be calculated within standard density functional 
theory and V is estimated from the projected density 
of states. The resonance energies as well as the non- 
adiabatic coupling parameters is obtained from excited 
state potential energy surfaces which are calculated with 
the method of linear expansion Delta Self-Consistent 
Field. The method is a generalization of standard Delta 
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FIG. 2: (Color online) The reaction yield of CO on Cu(lll) as 
a function of hot electron flux for four reaction energies. Note 
the initial linear dependence corresponding to single electron 
reactions. 
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System 


Experimental n 


Calculated n 


E R 


NO/Pd(lll) [2] 


3.3 


3.7 


1.0 


CO/Cu(lll) [4] 


3.7 


3.6 


0.4 eV 


CO/Cu(100) [7] 


8.0 ±1 


3.7 


0.5 


2 /Pt(lll) [5j 


3.0/5.6 


2.5/5.5 


0.4 e y 


0/Pt(lll) [8] 


15 


12 


0.8 eV 



TABLE I: Power law exponents obtained from (dJ-JS}. The 
systems involving NO, CO, and O2 are desorption experi- 
ments whereas 0/Pt(lll) is hot electron induced diffusion 
of atomic oxygen from a step to a terrace hollow site. For 
02/Pt(lll), the exponent depends on the laser frequency. 
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FIG. 3: (Color online) The yield as a function of electron flux 
obtained using the transition probabilities in equation (0 for 
different parameters (see text). 



code using the projector augmented wave method [21]. 
We modelled the surface by a three layer (4 x 4) super- 
cell with the top layer relaxed, a grid-spacing of 0.2 A, 
and a 4 x 4 surface K-point sampling. With the RPBE 
[22( | functional, we find Sq — 2.4 eV Hlu = 231 meV and 
A = —118 meV for the internal stretch vibration, and 
hw = 42 meV and A = —4 meV for the COM vibration. 
The internal mode completely dominates the transfer of 
energy from hot electrons to the molecule since it has 
a much larger coupling A, and the quantum of energy is 
five times larger than for the COM mode. Figure [1] shows 
the density of an excited top site molecule obtained with 
ASCF-DFT relative to the ground state density and one 
clearly sees the excited 7r* orbital. In figure [2] we have 
used (P)-© with hot electrons at e = 2.0 eV correspond- 
ing to the laser frequency used in [4| to calculate reac- 
tion rates, which require energies corresponding to 3, 4, 5 
and 6 internal vibrational excitations. In the non-linear 
regime corresponding to reactions induced by multiple 
scattering events, the rates are very well approximated 
by power laws with n = 2.8, n = 3.6, n = 4.3, and 
n = 5.1. Varying the parameters in the model reveals 
that in general, one always obtains very good power law 
fits with exponents n ~ ur. Increasing the energy of 
the hot electrons above the resonant energy eo tend to 
decrease n. 

Table U shows the calculated power law exponents for 
five systems compared with corresponding experiments. 
The DFT parameters use were the same as for CO on 
Cu(lll). We have assumed that a reaction occurs when 
an energy of 0.2 eV in excess of the reaction barrier has 
been transferred, which is consistent with measurements 
of the kinetic and internal energy distributions of des- 
orbed molecules [7[ and we have used hot electron ener- 
gies e = Sf + HQ where Q is the laser frequency. The 
agreement is very good except for CO on Cu(100). This 
could be due to the role of frustrated rotations, which 
has previously been found to couple strongly to metal- 
lic electrons in this particular system [23[. In general, 
however, the frustrated rotations have a low energy com- 
pared with the internal stretch mode and we expect their 
effect to be neglectable. Due to symmetry the frustrated 



rotations cannot couple linearly to the resonant electron 
and one would need a quadratic term like A2C a <4(6 + b^) 2 
in ([5j) to include these in the model [241 ]. 

A particularly interesting case is the hot electron in- 
duced desorption of O2 from Pt(lll) where a power law 
with exponent n ~ 5.6 ±0.7 was observed using a photon 
energy of 2.0 eV and an exponent ofn~3.0±0.5 using 
a photon energy of 4.0 eV [5]. The fact that the power 
law exponent depends on the photon energy contradicts 
the picture of a thermalized hot electron gas interacting 
with the molecule, which is the basic assumption in mod- 
els based on electronic friction. In contrast, the model 
(PQ)- ([5]) naturally gives rise to a decrease in the power law 
exponent when the energy of hot electrons is increased. 

The transition matrix P mn (e) has a very complicated 
structure and it is hard to extract the physics of the 
power law using these probabilities and (pQ)-(j4j). How- 
ever, the magnitude of P mn (£) is largely governed by the 
prefactor g n ~ m / (n — m)\ [17[ and in the following we will 
examine the consequences of assuming transition proba- 
bilities of the form 



Pmn — 6 



a' 1 



(n — m)V 



n > m, 



(6) 



where a is a dimensionless coupling constant, which de- 
scribes the coupling of hot electrons to the adsorbate vi- 
brational states. Repeated use of the algorithm (pQ)-(|4|) 
with these probabilities then reveals that to leading order 
in a one has 



a'" 



k-l 

3=0 



■jAt/T vib 



(7) 



We then consider a large flux e _A£ / T ^ b ~ 1 — At/T vib , 
sum up the geometric series, take the limit k — > 00 cor- 
responding to steady state, and get 



Q(n) = — (T m6 J ) n , 



(8) 



where At = 1/Jo- Thus, for small a the reaction proba- 
bility (|4]) will be dominated by such a term with n = ur. 

The power law emerges from summing up the detailed 
combinatorics of all possible ways of rising through the 
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vibrational states in the potential well. In figure [3^ we 
show the reaction yields for three values of ur and they 
are seen to approach power laws of the form Y oc Jq r 
for large fluxes. Even if a is not small, equations (pQ)-(|4j) 
tend to conserve the power law although the exponent 
becomes reduced from the value of ur when terms be- 
yond leading order are not vanishing. In figure [3)3 we 
show the yield when ur = 15 for a = 0.2, a = 0.5, 
and a = 1.0. At large fluxes the yields are well ap- 
proximated by power laws with exponents 14, 12, and 
10 respectively. One might worry that the fixed time 
interval At between scattering events is too crude an ap- 
proximation for the random nature of hot electrons inter- 
acting with the adsorbate. However, a sequence of time 
intervals with an average of At would lead to the 
replacement jAt — > YH=i tk-% m 0? which is well ap- 
proximated by jAt for large j. We have repeated the 
calculations leading to figure [3] but with the time inter- 
vals randomly drawn from an exponential distribution 
p(t) ~ e - £ / A ^ anc [ the power law is conserved on average. 

The interpretation of the power law exponent as the 
number of contributing vibrational states can be used to 
identify the reaction channel of a given adsorbate sys- 
tem. For example, in the study of hot electron mediated 
desorption of NO from Pd(lll) a power law with ex- 
ponent n ^ 3.3 was found. The internal stretch vibration 
corresponds to an energy of ftw ~ 210 meV whereas the 
other modes have vibrational energies ftw < 70 meV. 
Since the adsorption energy is E a ~ 1.0 eV we con- 
clude that the power law exponent ur ~ E a /hw has to 
arise from sequential excitation of the internal stretch vi- 
bration and subsequent anharmonic energy transfer to 
the desorption coordinate. In contrast, the study of 
hot electron induced surface diffusion of atomic oxygen 
on Pt(lll) Q gave rise to a power law with exponent 
n ^ 15, which fit very well with an experimental diffu- 
sion barrier of Ed — 0.8 eV and vibrational modes on the 
order huu ~ 50 meV. 

In summary, we have presented a theory of multiple 
inelastic scattering with transition probabilities calcu- 
lated in a non-adiabatic Newns- Anderson model, which 
lead to the ubiquitous power law of HEFatS and repro- 
duce experimentally found exponents. The interpreta- 
tion of the exponents as the number of contributing vi- 
brational states is a very useful tool to identify the re- 
action channel of a given system and also indicates that 
a classical treatmeant of the adsorbate motion is insuffi- 
cient. However, the model can only treat quadratic po- 
tentials and are thus not able to treat the anharmonic 
effects, which are expected to play an important role in 
the transfer of internal vibrational energy to the reac- 
tion coordinate. In the case of molecular desorption the 
model applies because the internal stretch mode, which 
is well approximated by a harmonic potential, dominates 



the energy transfer. Furthermore, while our assump- 
tion of constant hot electron energy e = Sf + HQ gives 
the right dependence of the exponents on Q, the hot 
electrons proceeding a laser pulse will undergo scatter- 
ing and produce some distribution of electron energies. 
The true (time-dependent) distribution lies somewhere 
between the present assumption and a thermalized hot 
electron gas and the model should thus be regarded as 
complementary to a statistical approach based on elec- 
tronic friction and Langevin dynamics, which assumes a 
thermalized hot electron gas. 
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